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ABSTRACT 


This thesis is devoted to the study of standard perturbation 
techniques applied to two well known second order differential operators. 
In the introductory chapter we present the definitions and basic theorems 
that are necessary to develop these techniques. Chapter two deals with 
Legendre's equation which is perturbed with two unbounded operators. 

In the third chapter we study Bessel's transformed equation perturbed 

in the standard way. One example shows the validity of these techniques 
for any bounded perturbing operator and the second example is a differ- 

ential perturbing operator. The last chapter deals with some perturbing 


operators that produced no interesting results. 
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CHAPTER I 


INTRODUCTION 


Perturbation techniques have been used on ordinary differential 
euqations by Bellman [1], Rellich [1], and Cole [2]. In this thesis we 
use perturbation techniques to look for series solutions to perturbed 
forms of Legendre's equation and Bessel's equation. Some particular 
examples of a perturbing operator are investigated. In the last chapter, 


we indicate some cases for which no results were obtained. 


1. PRELIMINARY REMARKS. 


The vector space of real valued functions that are continuous 
for all xe [a,b] is denoted by [a,b] and the space of real valued 
functions that are twice continuously differentiable for all x « [a,b] 
is denoted by Geieen) . Finally, real valued functions that are 
continuous at all but a finite number of points of an interval and at 
each point of discontinuity, both the left and right limits of the 
functions exist, are called piecewise continuous real valued functions. 


We denote this vector space by P@[a,b] 


Let the real vector space V _ be defined such that for any 
f;2 € V, the inner product 
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(f,2) = i £ (x) g(x) dx Cie.) 
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6éxists for x P [a,b] . V is called a Euclidean or inner product space. 
The vector space @* [ash] with this inner product is denoted by Cte 
The space P@[a.b], again with (1.1) as inner product, is the Euclidean 


space PC[a,b] and PC[a,b] contains C*Ta,b] as a subspace. 


For any fe¢«V, the non-negative real number 


lel] = V4) (189) 
is the norm of f. A sequence fi of5s--- of real valued functions 


of V is said to be orthogonal if Cie pe) = 0 whenever i + j and 
[£5 ll 0 for every j . It is orthonormal if it is orthogonal and 
Ile, ll 1 for every j . As we shall see later, Legendre polynomials 


and Bessel functions form orthogonal sequences which are easily normalized. 


Let the sequence if} for f, Gat Wore 2 stress Dew all 


orthonormal sequence of functions. Also let f¢«V.. The series 


fl 148 


froin (133) 
; 1 ali 

i=1 
is said to be the orthogonal development of f in V with respect to 
the f,'s . We say that the sequence if} wl Sees care (Complete 
if for all f ¢«V , the orthogonal development of f converges to f 
This will be shown to be the case for sequences of Legendre polynomials 


and Bessel functions in PC[-1,1] and PC[0,1] respectively [10]. 


We are concerned with Legendre's differential equation 
(Chapter II) and the transformed Bessel's differential equation 


(Chapter III). In each case we have an equation of the form 


(S58) 
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Ly es eAy (1.4) 
where the differential operator L is defined by 


Ly = (p(x)y')’ + aq(x)y , xe [a,b] . (iS) 


An operator of the form (1.5) is said to be self-adjoint. The 
necessary and sufficient condition for the general second order differ- 


ential operator 
Dy = fy'' + gy' + hy (1 26) 


Cinawhich ieee yeand i) ssarerfunctions of. x)» to be self-adjoint is that 
e= f°. [ll]. This ds clearly the case for Leeéendre's differential 


operator where g = f' = -2x and for Bessel's differential operator where 


The values of i which yield nontrivial solutions of (1.4) 
are called eigenvalues and the assoicated solutions are called eigen- 
functions or eigenvectors. If L is a linear self-adjoint operator then 


all of the eigenvalues are real[10]. 


The set of eigenvalues is called the spectrum [7]. Titchmarsh 
[LG Chi... -3)| Soiree some criteria that enable us to decide when the problem 
(1.4) has a discrete or continuous spectrum. The most general conditions 
for discreteness or continuity are unknown. In this thesis we assume, 
as is the case for Legendre's equation and Bessel's equation, a discrete 


spectrum which is bounded above. 


Dunford and Schwarz [6, p. 1478] prove that a second order self- 
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adjoint operator with a discrete spectrum that is bounded above or below 
has non-degenerate eigenvalues. Therefore, we can disregard the problem 


of degeneracy. 


2 
Let S be a subspace of C [a,b] determined by a pair of 


boundary conditions such that 


|b 
1 ae 1 = 
p(x) Ly, (x)y, Ge) y, (xy (x)] : 0 
for every xy, € SS. Then any set of eigenvectors belonging to distinct 
i 2 


eigenvalues for the problem 
Deano. 
is orthogonal in C[a,b] with the inner product 


b 
(V4 Vo) ee ha ¥1 C&) ¥ 5 (x) dx 


Also 
(Ly, .¥5) = (y, Ly.) 


for Ya: Y>5 vectors of a Euclidean space [10]. 


2. THE PERTURBATION PROBLEM. 


The method of perturbation [1], [2], [7] is applicable if the 
operator to be analyzed is, in one sense or another, near to an operator 
whose spectrum is known. The spectrum of an operator is sensitive even 


to the slightest changes in the operator. 
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The operator with a known spectrum, the "unperturbed" operator, 


will be denoted by L , the "perturbing" operator by B , and the 


"perturbed" operator depending on a parameter ce is denoted by L + eB 


We assume, as is the case for Legendre's differential operator 
and Bessel's differential operator, that the unperturbed operator L , 
given by (1.5), is self-adjoint and has discrete isolated eigenvalues 


doi with corresponding normalized eigenvectors vee 


Suppose the perturbed operator L+ eB has an eigenvalue M 
which, together with an eigenvector veer depends analytically on e 


Then expansions 
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2 
A. = Yo + Ye + Yo€ ore (1.9) 


may be valid for sufficiently small values of ¢ and may be formally 


substituted into 


(LteB)y = Rian : C1210) 


We equate coefficients of like powers of e . In the case that the 
eigenvalues are distinct, if other orders of « are used, the correspond- 
ing Y 5 would be identically zero. The following sequence of equations 


is obtained: 


(L-r) Yq =.0 , 
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Of course, the vector Vie is not uniquely determined by these 
equations since it may be multiplied by any function of e« . Since 
Yo + 0 this factor could be so chosen that the norm lly. | of ye is 


1 for sufficiently small e¢ [15]. Instead, we impose the condition 
which implies [7] 
(Yo°¥4) = (Yo V5) im Céokicn oS 0 ° (lee 2) 


Now we can establish the following theorem. 


THEOREM 1.1. Assuming a point spectrum of non-degenerate eigenvalues, 


the right member of each equation in (1.11) is orthogonal to Yo: 


PROOF: We take the inner product of the unperturbed equation with Ys > 


i= Deets) BOLNCe L-r) is a linear self-adjoint operator, 
((L-r> 0) 95 2¥4) = (yoo (L-Aj)y5) =O, ce 


Therefore the left hand side of the system (1.11) is orthogonal to ae 


which yields the desired result. 


De do is an isolated eigenvalue of L and the domain of 
L-) is the space spanned by the ortho-complement set YpoYoores of 
Yo then (L-r > exists. Therefore the sequence of equations (1.11) 
can be solved for Yoo¥po°°* (the first being satisfied by definition) 
and Yi2Yoo+++ are made unique by (1.12). Since the sequences ito : 


n= 0,1,..., that we deal with in Chapters II and III are complete [10] 
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Theorem (1.1) and the equality (1.12) yield the values 
: A = (BY, »VQ) > é do a (BY) Yo) o+-> ele) 
of the expansion coefficients of the eigenvalue d. 


A linear operator L , which assigns to each vector y in 
V avector Ly in V, is called "bounded" if the inequality 
|Ly| < aly| holds for all vectors y with an appropriate number a. 
If B is a linear self-adjoint bounded operator in a neighbourhood 
of e€ = 0 , the convergence of (1.8) and (1.9) is guaranteed by 
Rellich [12, p. 57]. We study a bounded operator in Chapter III. How- 
ever, if this is not the case, although the formal power series may 
always be written down, the convergence of these series must be 
demonstrated explicitly in each case, as we will do with examples in 


Chapters ILI and III. 
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CHAPTER II 


LEGENDRE'S DIFFERENTIAL EQUATION 


In this chapter, we study Legendre's differential equation 
[14, p. 76] using the perturbation techniques of Chapter I. Some 
intersting results are obtained, particularly when specific examples 
of the perturbing operator B are used. When a simple first order 
differential operator is chosen, it is found that only the eigenvectors 
have nonzero perturbing terms. However, when a situation similar to 
Legendre's associated equation [14, p. 297] is examined, we find that 
both the eigenvectors and the eigenvalues have nonzero perturbing 


terms. 


1. LEGENDRE POLYNOMIALS. 


The Legendre polynomials may be defined by Rodrigues' formula 
[ee p+ P201)% 
n 


=a . (x°-1 
2 ni <ax 


P(x) = Pee sae OC ok, (oui) 


It is clear that for each fixed n , P fs) is a polynomial of degree 
no; P fs) contains only odd powers of x when n is odd and only 


even powers when n is even. 


Lemma 2.1. Pf) ®akyx P f-)) = (one for all on. 
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PROOF: The equality PD) = 1 obviously holds for P) 1° 


Moreover, if we assume that Py (LD) a PS) oy 20.5 Lesh seche 


recurrence relation [10, p. 416] 
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ser ( Tah ni ea ei (2.2) 
gives 
2n+1 n 


The desired result follows by mathematical induction. In much the same 


way it can be shown that Pf) Se anya a 


We frequently need to compute integrals of the form 


1 
10 fC) Pi (x).dx 
Lap 


If f and its first n derivatives are continuous throughout the 


interval =1 <x <1, we use integration by parts n times. Since 


: . 2 : : 
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THEOREM 2.1. If £ and its first n derivatives are continuous through- 
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dl 
f sGop (dx = £2 fy ®— swax . (2.3) 
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F ‘ m f . 
In particular, 1f£° f(x)"= x ,* where’ ~m is an integer 


OF amt<: 1 


» we have 


1 
fo xP Gx)dx= 0, Mesto nal. (2.4) 
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2. LEGENDRE'S DIFFERENTIAL EQUATION. 


The boundary conditions for Legendre's differential equation 


oo . oe moe SU os Shs Cela) 


where Lys ((1-x")y!)' » are simply: the solutions remain finite at 
x =+1. When ro = -n(ntl) , n= 0,1,2,..., the Legendre polynomials, 


P (x) , are the only bounded solutions of (2.5). 


It is easy to prove [3, p. 202] that the sequence of functions 


{P(x} > m= 0,122....,. ts )an-orthogonal sequence for, x « [-1,1]. 


THEOREM 2.2. For non-negative integral values of n and m, xe [-1,1] 


and the inner product 
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CaN Lge P (x)P_(x)dx , (2.6) 
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for an arbitrary function f in PC[-1,1] , is known as the Legendre 
series expansion of the function f . The proof of the following theorem 


is on page 423 of Kreider [10]. 
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Clearly, each successive eigenvalue will contain an expression 


of the form: 
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PROOF: Using integration by parts, the case when i =n follows 
immediately. If i> nn, the result follows from Theorem 2.1. If 
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CHAPTER III 


BESSEL'S DIFFERENTIAL EQUATION 


Chapter III deals with Bessel's differential equation [10]. 
Once we have written the equation in self-adjoint form we use standard 
perturbation techniques [7] to obtain a system of equation [10]. Two 
examples, one with a general bounded linear perturbing operator and one 


with an unbounded linear perturbing operator, are studied. 


1. THE ZEROS OF SOLUTIONS OF BESSEL'S EQUATION. 


In physical applications Bessel's equation usually arises in 


the self-adjoint form 
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only impose the boundary condition y(1) =0. 


The following theorem is proved [10, p. 235] by applying 


Sturm's comparison theorem [10, p. 232] to 
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which is obtained from (3.1) by making the change in variable y = z/vx ° 
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and the theorem is proved. We shall denote the right side of this 


inequality by K(m) 
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are eigenvectors. Since x does not vanish when x=1, k =0 


is not an eigenvalue. 
Since we will soon be computing series expansions relative 
to the eigenvector found above, we now evaluate their norms in PC[0,1] 
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Now that we have stablished the existence of an infinite set 
of orthornormal functions for the boundary value problem (3.4) we need 
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2. PERTURBATION OF BESSEL'S EQUATION BY A BOUNDED OPERATOR. 
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CHAPTER IV 


CONCLUSION 


Perturbation of eigenvalue problems has been thoroughly 
studied, particularly for bounded perturbing operators, by Franz Rellich 
[12], [13]. His fundamental results [13] have not been generalized 
for the unbounded case. The application of these techniques to 
Legendre's second order differential operator [10] and Bessel's second 


order differential operator [10] has not been found in the literature. 


When we finally set these problems up correctly, we were left 
with the difficult task of choosing perturbing operators that would 
yield meaningful results. Rellich shows that a small parameter e 


does not necessarily imply a small perturbation [12]. 


With respect to Legendre's differential equation, we 


began by considering the equation 
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anticipated, relatively uninteresting results were obtained - the 
eigenvectors were not perturbed but the eigenvalues were. Finally 
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finite sums and could be easily bounded. 


Several attempts at perturbing Bessel's equation were 
unsuccessful. We considered 
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expressions for this problem contained x and were not explicitly 

integrable. We were unable to find any estimates on regions of 


convergence for this case. 


The transformation y = z/Vx gave us a nice boundary value 
problem to perturb. First of all, we considered a bounded per eateine 
operator to demonstrate Rellich's result that a convergent series 
exists for the eigenvalues and eigenvectors of the perturbed problem. 
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In this thesis, we have chosen the perturbing operators B 
with no obvious physical importance. However, the fact that first 
order differential perturbing operators yield well defined eigenvalues 


and eigenvectors for the perturbed eigenvalue problem is significant. 
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